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REMOTE  INFRARED  VAPOR  DETECTION  OF  VOLATILE  ORGANIC  COMPOUNDS 


1 .  INTRODUCTION 

Environmental  monitoring  with  Fourier  transform 
infrared  (FT-IR)  spectrometry  offers  a  robust  method  for  the 
identification  and  quantification  of  various  gas  pollutants. 
Volatile  organic  compounds  (VOCs)  represent  a  major  pollutant 
from  industrial*,  automotive^,  and  waste  sight’  sources.  These  VOC 
sources  have  been  characterized  by  open  path  air  monitoring  with 
a  variety  of  FT-IR  spectrometer  configurations.  Two  spectrometer 
configurations  are  active  bistatic  and  passive  (see  Figure  1) 

The  active  bistatic  configuration  uses  an  infrared  (IR) 
source  and  provides  transmission  through  the  target  vapor  to  the 
FT-IR  spectrometer.  This  configuration  resembles  the  conventional 
laboratory  spectrometer  with  two  possible  exceptions.  First,  the 
sample  cell  is  replaced  with  an  open  path  approximately  100 
meters  in  length.  Unlike  the  controlled  laboratory  cell,  the  gas 
sample  parameters  are  determined  by  meteorological  conditions.^ 
Second,  the  sample  is  located  between  the  IR  source  and 
interferometer,  rather  than  between  interferometer  and  detector. 
Thus,  operation  of  the  FT-IR  spectrometer  in  the  passive 
configuration  only  requires  substitution  of  a  terrestrial  scene 
for  the  IR  source. 

The  passive  FT-IR  spectrometer  configuration  relies  on 
a  radiance  difference  between  the  target  vapor  cloud  and 
background  for  an  emission  spectrum.  The  radiance  is  the  spectral 
radiant  emittance  scaled  by  an  emissivity  factor.’  The  spectral 
emittance  is  calculated  from  Plank's  function  and  depends  on 
temperature.  Thus,  the  radiance  difference  can  be  interpreted 
primarily  as  a  temperature  difference  between  target  vapor  cloud 
and  background.*  Possible  background  scenes  are  terrestrial,  sky, 
or  some  mixture  of  both  terrestrial  and  sky.  These  background 
scenes  exhibit  large  temporal  and  spatial  variations.  The 
variations  are  caused  by  the  spectral  dependence  of  emissivity 
values  and  temperature  fluctuations  in  the  atmosphere.  The 
background  variations  preclude  the  collection  of  a  stable 
background  reference  spectrum. 

This  investigation  considers  signal  processing 
approaches,  which  do  not  rely  on  the  existence  of  a  stable 
reference  background.’***  The  approach  suppresses  the  broad  band 
detector  envelope,  while  maintaining  the  narrower  band  spectral 
signatures  of  the  target  vapor  cloud.  The  broad  band  suppression 
requires  three  steps.  First,  selection  is  made  of  an 
interferogram  segment  where  the  narrow  band  target  signature 
remains  relatively  strong  compared  to  the  attenuated  broad  band 
background  signal.**  Second,  a  bandpass  digital  filter  is  applied 
to  the  selected  interferogram  segment  to  remove  signals  at 
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(a) 


(b) 


(0 


Figure  1.  Experimental  Methods  for  Data  Collection: 

(a)  Active  Bistatic,  (b)  Passive  Terrestrial,  and 
(c)  Passive  Laboratory  Spectrometer  Configurations. 


Figure  2.  Flow  Chart  of  the  Data  Collection  Program. 
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frequencies  other  than  those  corresponding  to  a  selected  band  of 
the  target  vapor.  Third,  pattern  discrimination  techniques  are 
applied  to  the  filtered  interferogram  segment.  These 
discrimination  techniques  effectively  examine  the  narrow  band 
shape  of  the  target  vapor  to  avoid  erroneous  detections  due  to 
interferant  vapors.  The  three  step  background  suppression  process 
is  demonstrated  for  two  VOCs  (i.e.,  acetone  and  methyl  ethyl 
ketone)  as  well  as  a  narrower  band  test  gas  (i.e.,  sulfur 
hexafluoride) .  Interferogram  data  is  collected  under  a  variety  of 
experimental  conditions  to  evaluate  the  effectiveness  of  the 
three  step  background  suppression  technique. 


2 .  INSTRUMENTATION 

The  MIDAC  FT-IR  emission  spectrometer  (serial  number 
120)  with  a  Newtonian  telescope  and  remote  IR  source  is  used  in 
this  investigation  (Midac,  Inc.,  Irvine,  CA) .  The  FT-IR 
spectrometer  consists  of  a  linear  drive  Michelson  interferometer 
with  a  germanium  coated  zinc  sr lenide  beamsplitter/compensator 
plate.  An  off-axis  parabolic  mirror  focuses  the  IR  radiation  onto 
a  2x2-mm  liquid  nitrogen  cooled  Hg:Cd:Te  detector  that  is 
designed  for  8  to  12  micron  response.  The  spectrometer  is  found 
to  possess  a  6°  field  of  view  (FOV) .*  The  Newtonian  telescope 
provides  a  limited  0.3®  (FOV)  with  a  500  cm^  aperture  which  can 
be  filled  with  the  IR  source.  The  IR  source  is  composed  of  a 
2000  cm^  aperture  metal  paraboloid  to  collimate  a  1550  K  silicon 
carbide  heater  element  point  source.  The  interferometer  scan 
speed  is  established  to  be  2.45410.005  mm/sec  for  these 
measurements.'^'^  The  interferogram  data  is  sampled  at  every 
eighth  zero  crossing  for  a  total  of  1024  sampled  points.  Spectra 
that  are  computed  from  these  interferograms  have  a  spectral  point 
spacing  of  approximately  4  cm'. 

The  spectrometer  is  interfaced  to  a  Dell  System  486P/50 
IBM  PC  compatible  computer  operating  under  MS-DOS  Version  5.0 
(Microsoft,  Redmond,  WA) .  Interferogram  collection  is  performed 
with  the  MIDCOL  software  package.'*  This  software  is  written  in 
Microsoft  C,  Version  6.1  and  utilizes  graphics  routines  from 
Essential  Graphics  System  (South  Mountain  Software,  Inc.,  Orange, 
NJ) .  The  program  permits  rapid  data  collection  of  interferograms 
in  16-bit  integer  format  as  delivered  from  the  spectrometer's 
analog-to-digital  converter.  Figure  2  provides  a  flow  chart  of 
the  data  acquisition  portion  of  the  MIDCOL  software  and  outlines 
its  functions.  The  MIDCOL  software  places  interferograms  into 


*Tveten,  L.H.,  "Acceptance  Test  Procedure  for  the  Background 
Measurement  Spectroradiometer , "  Technical  Report  CDRL  Sea.  No. 
A019 .  ATP  25368-01,  DAAK-79-C-0051 ,  Honeywell  Inc.,  St.  Petersburg, 
FL,  October  1979. 
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(a)  BACKGROUND 


<a)  BACKGROUND 


Figure  3 . 


Figure  4. 


Single-Beam  Passive  Terres¬ 
trial  Spectra  with  Comparable 
Responses:  (a)  Background, 

(b)  Acetone,  and  (c)  Acetone/ 
Background  Difference. 


Single-Beam  Passive  Terres¬ 
trial  Spectra  with  Dissimilar 
Responses:  (a)  Background, 

(b)  Acetone,  and  (c)  Acetone/ 
Background  Difference. 
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aggregate  data  files  that  are  composed  of  up  to  3000 
interferogram  subfiles.  The  data  files  are  prefixed  with  a  global 
header  which  contains  filename,  date,  time  of  day,  interferometer 
configuration,  meteorological  data,  comments,  and  number  of 
interferogram  subfiles.  Each  interferogram  subfile  is  identified 
by  a  subfile  header  which  includes  scan  number,  time  of  day, 
center  burst  location,  and  error  code.  This  header  information 
allows  a  systematic  labeling  and  categorization  of  the  data 
important  in  subsequent  data  analysis.  Appendix  A  contains  the 
file  format  specifications. 


3 .  EXPERIMENTAL  METHODS 

The  limit  of  detection  for  remote  FT-IR  spectrometers 
is  governed  by  four  parameters:  (1)  radiance  difference  between 
IR  background  source  and  target  vapor,  (2)  vapor  path  length,  (3) 
vapor  concentration,  and  (4)  atmospheric  attenuation.  These 
parameters  are  present  for  all  three  experimental  methods  that 
are  used  in  this  investigation  (see  Figure  1) .  The  active 
bistatic  configuration,  shown  in  Figure  la,  provides  an  IR  source 
of  such  intensity  that  the  lower  vapor  radiance  becomes 
insignificant.  A  transmission  experiment  is  performed  when  the  IR 
source  is  operated  at  1550  K.  Distances  from  the  IR  source  to  the 
FT-IR  spectrometer  varied  from  25-100  meters.  In  the  bistatic 
configuration,  variations  of  the  vapor  cloud  concentration  and 
atmospheric  attenuation  control  the  limit  of  detection.*’  If  the 
IR  source  temperature  is  significantly  reduced  and/or  the  IR 
source  does  not  fill  the  spectrometer  FOV,  then  the  IR  source 
radiance  becomes  comparable  to  that  of  the  vapor.  Under  these 
conditions,  the  spectral  response  of  the  bistatic  configuration 
resembles  that  of  the  passive  approach. 

The  passive  terrestrial  approach  is  described  in 
Figure  lb  and  uses  only  the  ambient  background  terrain  radiance. 
The  1550  K  IR  source  provides  an  initial  alignment  and  is 
subsequently  removed  from  the  spectrometer  FOV.  In  the  passive 
case,  the  radiance  difference  between  the  ambient  background 
(i.e.,  about  300  K)  and  the  target  vapor  is  often  small.  This 
small  radiance  difference  is  usually  the  controlling  factor  in 
detection.  If  the  contrast  radiance  between  the  vapor  cloud  and 
terrain  background  is  insufficient,  then  no  spectral  target  vapor 
detection  is  possible. For  the  bistatic  and  passive 
configurations,  VOC  vapors  have  been  dispensed  by  either  static 
or  forced  air  flow  evaporation  of  spectral  grade  Aldrich 
chemicals.  Sulfur  hexafluoride  (SF^)  releases  have  been  made  with 
a  fixed  flow  rate  from  a  Matheson  Gas  Products  QA  gas  cylinder.'* 

The  third  and  final  approach  in  Figure  Ic 
experimentally  simulates  the  passive  terrestrial  measurements  in 
the  laboratory.  The  telescope  is  not  used  in  this  experiment.  The 
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vapor  source  is  a  small  aperture  and  short  path  length  vapor 
cell.'^  The  short  path  vapor  cell  provides  an  aperture  of  62  cm^ 
and  a  path  length  of  8.2  cm.  The  4x4-inch  extended  SR-80 
blackbody"  furnishes  an  adjustable  temperature  source  from  5  to 
50  *C.  The  blackbody  is  a  NIST  traceable  source  that  is  accurate 
to  ±0.03  *C  and  precise  to  ±0.01  •€.  The  optical  train  of  FT-IR 
spectrometer,  vapor  cell,  and  blackbody  is  aligned  with  a  helium 
neon  laser.  Distances  between  spectrometer /cell  and 
cell /blackbody  are  approximately  10  cm.  These  distances  along  the 
optical  train  axis  ensure  that  the  spectrometer  FOV  lies  within 
the  62  cm^  aperture  of  the  short  path  cell  and  views  only  the 
4x4-inch  blackbody  IR  background. 

A  variety  of  conditions  are  used  in  data  collection 
with  the  spectrometer  configurations  of  active  bistatic,  passive 
terrestrial,  and  passive  laboratory  that  are  illustrated  in 
Figure  1.  Table  1  summarizes  the  number  of  interferograms 
collected  with  each  configuration  for  each  target  vapor.  Separate 
training  and  prediction  data  sets  of  interferograms  have  been 
acquired  for  each  vapor  (see  Table  1) .  The  training  set 
interferograms  are  used  in  selection  of  the  appropriate 
interferogram  segment,  computation  of  the  optimal  bandpass 
digital  filter,  and  optimization  of  the  pattern  discrimination 
procedures.  The  prediction  set  interferograms  serve  as  an 
independent  test  set  for  evaluation  of  the  background  suppression 
procedure.  The  complete  training  set  contains  23510 
interferograms  for  acetone,  22960  interferograms  for  MEK,  and 
25414  interferograms  for  SF^.  The  complete  prediction  set  holds 
45978  interferograms  for  acetone,  48804  interferograms  for  MEK, 
and  27124  interferograms  for  SF^.  The  representative  background 
interferograms  for  the  training  and  prediction  subsets  of  each 
vapor  are  determined  by  the  global  sampling  algorithm  developed 
by  Carpenter  and  Small.’’  Numbers  of  interferograms,  that  are 
selected,  are  listed  in  Table  1.  Approximately  1000 
interferograms  for  the  training  subsets  of  each  compound  are 
selected  on  the  basis  of  visual  inspection  of  the  Fourier 
transformed  spectra  (see  Appendix  B  for  types  of  interferograms 
selected) .  Similarly,  a  number  of  vapor  containing  interferograms 
are  selected  for  the  prediction  subsets  giving  a  total  of  10000 
interferograms  in  each  subset.  The  passive  laboratory 
interferograms  have  provided  an  accurate  test  to  verify  the 
visual  inspection  of  spectra  for  conditions  near  the 
spectrometer's  detection  limit. 


4 .  ANALYSIS 

4.1  Fr^q^e.ocY-B.Qffiain  Analysis. 

Spectral  analysis  for  vapors  in  the  IR  depends  on  the 
parameters  of  band  position,  width,  strength  and  contour.  Table  2 
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summarizes  both  the  spectra  1“**'*“  and  vapor  pressure^  information 
pertinent  to  the  vapors  that  have  been  investigated  in  this 
study.  Columns  three  and  four  of  Table  2  contain  the 
absorptivities  and  vapor  pressures.  These  provide  an  estimate  of 
relative  band  intensities  assuming  an  uniform  detector  response. 
Unfortunately,  photoconductor  detectors  usually  possess  an 
uniform  response  over  very  limited  spectral  ranges  as 
demonstrated  in  Figure  3a  for  a  single-beam  passive  background 
spectrum.  The  band  positions  for  the  VOCs  and  SF^  are  also 
delineated  in  the  background  spectrum  of  Figure  3a  and  tabulated 
in  column  one  of  Table  2.  The  band  widths  are  listed  in  column 
two  of  Table  2  as  the  full-width-at-half-height  (FWHH)  values.  A 
band  overlap  of  about  100  cm*  occurs  between  acetone  and  methyl 
ethyl  ketone  (MEK) .  In  addition,  the  detector  response  for  these 
VOC  band  positions  is  only  about  one  half  that  of  the  SF^  band 
position. 


The  lower  detector  response  combined  with  a  decrease  in 
the  VOC  absorptivity  by  a  factor  of  20  to  40  as  compared  to  SFg 
makes  the  detection  of  VOCs  much  more  difficult  than  SF^.  The 
detection  is  further  complicated  by  variation  of  the  background 
for  passive  terrestrial  measurements.  Figures  3  and  4  show  the 
difference  between  an  acetone  spectrum  and  two  separate 
background  spectra.  The  background  and  acetone  spectra  in 
Figure  3  are  collected  under  similar  conditions.  The  difference 
spectrum  contains  a  prominent  acetone  feature.  Measurement  of 
another  background  only  a  few  minutes  later  shows  a  substantial 
change  in  the  background  radiance  has  occurred.  In  this 
difference  spectrum  (see  Figure  4c) ,  the  acetone  feature  becomes 
hidden  by  a  large  baseline  fluctuation.  This  emphasizes  the 
potential  problem  with  attempts  of  averaging  consecutive  spectra 
to  increase  the  signal-to-noise  ratio  (SNR) .  For  spectra  with  low 
SNRs,  the  background  variation  also  limits  the  utility  of  the 
traditional  approach  of  a  reference  background  subtraction 
followed  by  classification  with  a  linear  discriminant.® 
Unfortunately,  it  is  often  impossible  to  obtain  an  accurate 
reference  background  spectrum.  Therefore,  an  alternate  method  of 
background  removal  or  suppression  is  necessary  for  this 
application. 

4.2  Time-Domain  Analysis. 

The  superposition  of  a  set  of  damped  cosine  waves 
describes  an  IR  inter ferogram.  The  set  of  damped  cosine  waves  is 
defined  by  the  IR  radiation  detected  with  the  spectrometer.  The 
cosine  damping  for  each  frequency  defines  the  composition  of 
frequencies  distributed  across  the  inter ferogram."  A  large 
damping  factor  indicates  a  wide  band  in  the  frequency  domain, 
while  a  small  damping  factor  corresponds  to  a  narrow  band.  Both 
the  digital  filtering  and  the  Gram-Schmidt  vector 
orthogonal ization^  procedures  for  direct  interferogram  analysis 
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rely  on  the  broad  band  of  the  detector  response  damping  out 
before  the  narrow  band  of  the  target  compound.  This  effect 
results  in  the  selection  of  interferogram  segments  that  are 
displaced  from  the  zero  path  difference  (ZPD)  of  the 
interferogram  center  burst. A  comparison  of  the  digital 
filter  and  Gram-Schmidt  approaches  by  Bjerga  and  Small"  indicates 
that  the  digital  filter  achieves  a  lower  SNR  detection  and 
provides  enhanced  computational  efficiency.  Therefore,  the 
digital  filter  approach  is  used  for  the  remote  monitoring  of 
acetone  and  MEK  which  have  SNRs  significantly  lower  than  SF^.  The 
generation  of  the  digital  filter  requires  the  judicious  selection 
of  both  the  filter  type  and  interferogram  segment.  Subsequent 
examination  of  the  digitally  filtered  interferogram  segments  with 
pattern  recognition  techniques  furnishes  a  discrimination 
capability  which  depends  on  the  spectral  band  contour  of  the 
target  vapor. 


Table  2.  Spectral  and  Vapor  Pressure  Information 


^peak 

(cm') 

FWHH 

(cm') 

a(lO') 

(mVmg) 

p.^@2rc 

(mm  of  Hg) 

Acetone 

1217”” 

1216+ 

49+ 

1.3” 

2 . 0” 

193.9” 

MEK 

1175” 

1172+ 

29+ 

0.5” 

74.5” 

SFs 

939” 

947” 

946+ 

10+ 

47.9” 

t  From  NIDAC  spectrometer  operated  at  4  cm*'  resolution. 


4.2.1  Digital  Filtering. 

Digital  filtering  is  a  standard  means  of  removing  unwanted 
noise  frequency  components  by  smoothing.^  Digital  filters  also 
provide  the  ability  to  remove  other  selected  interference 
frequencies.”  The  most  common  type  of  digital  filter  for 
selective  frequency  attenuation  is  the  finite-impulse-response 
(FIR)  filter.  Equation  1  is  the  functional  form  of  the  FIR 
filter. 

Yi*  =  +....+  fN+i»Yi  (1) 

Yj*  denotes  the  resultant  filtered  data  point  •  Yi.N  to  Yi  represent 
the  input  unfiltered  data  points.  The  coefficient  weights,  to 
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f,4^,  determine  the  filter's  frequency  response.  If  each 
coefficient  weight  is  set  to  1/(N+1),  then  a  moving  average 
results.  The  moving  average  filter  has  a  low  pass  frequency 
response.  Digital  filtering  for  the  extraction  of  specific 
spectral  bands  in  Interferograms  requires  a  response  that  depends 
on  both  the  spectral  band  and  detector  response  curve. 

A  FIR  filter  has  previously  been  developed  to  extract 
the  SF«  spectral  feature  at  940  cm'  from  interferograms  collected 
with  a  passive  FT-IR  spectrometer  mounted  on  a  helicopter The 
derivation  of  the  filter  coefficients  begins  in  the  frequency 
domain.  First,  a  multiplication  of  a  Gaussian  function  with  the 
spectrum  containing  SF«  is  performed.  The  Gaussian  function  is 
centered  on  the  SF«  band  at  940  cm'  and  possesses  a  bandwidth  of 
14  cm'.  The  multiplication  results  in  an  ideal  filtered  spectrum. 
The  inverse  Fourier  transform  provides  the  filtered 
interferogram.  To  perform  the  filtering  function  in  the  time 
domain,  this  Investigation  uses  equation  1.  Equation  l  is  a 
linear  model  that  relates  the  set  of  independent  variables  (raw 
interferogram  points)  and  dependent  variable  (filtered 
interferogram  points) .  The  coefficient  weights  are  determined  by 
regression  analysis  on  a  set  of  interferograms  containing  SF^. 

The  constraint  is  imposed  on  the  regression  analysis  that 
consecutive  interferogram  points  must  be  used  in  generation  of 
each  filtered  interferogram  point. 

Subsequent  work  has  focused  on  removal  of  the 
constraint  requiring  that  each  filtered  interferogram  point  use 
the  same  number  of  consecutive  raw  interferogram  points."  This 
permits  a  different  filter  model  to  be  used  for  each  point.  A 
stepwise  multiple  regression  permits  identification  and  deletion 
of  the  filter  coefficient  terms  that  are  less  significant  in 
approximation  of  the  ideal  filtered  interferogram.  This  type  of 
filter  is  called  a  finite  impulse  response  matrix  (FIRM)  filter. 
The  fewer  filter  coefficient  terms  of  the  FIRM  filter  for  a  given 
interferogram  point  reduces  the  dynamic  range  requirements  that 
are  placed  on  the  filter  coefficients.  This  has  given  an  improved 
performance  in  computational  efficiency  and  SNR  enhancement. 

This  investigation  used  FIRM  filters.  The  optimal 
interferogram  segment  location  for  the  narrow  band  SFj  was  a  75 
point  segment  that  was  170  interferogram  points  displaced  from 
the  ZPD.  A  FIRM  filter  was  generated  using  this  segment  for  the 
training  set  of  928  interferograms  containing  SFj  and  9072 
background  interferograms.  The  percentage  distribution  of 
interferograms  from  the  active  bistatic,  passive  terrestrial,  and 
passive  laboratory  configurations  was  24:59:17  respectively  (see 
Table  1  coliunns  3  and  4) .  Coefficient  terms  were  selected  from 
points  i-0  to  i-100  in  equation  1  using  stepwise  regression.  The 
entrance  criterion  for  coefficient  term  inclusion  in  the  filter 
was  at  a  significance  level  of  99.99%  based  on  the  F  test.  The 
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average  value  for  the  regression  calculations  was  99.0%  and 
the  average  number  of  coefficients  required  to  filter  each  point 
was  thirty. 

The  application  of  the  FIRM  filter  is  illustrated  in 
Figure  5.  The  Figures  5a  and  5b  are  plots  of  raw  interferogram 
segments,  while  Figures  5c  and  5d  are  plots  of  the  corresponding 
FIRM  filtered  segments.  Figures  5a  and  5c  show  the  effect  of  the 
filter  on  a  background  interferogram  segment  with  SF^  absent. 
Figures  5b  and  5d  demonstrate  the  filtered  response  when  SF«  is 
present  in  the  raw  interferogram  segment.  Clearly,  Figure  5 
illustrates  the  effectiveness  of  the  FIRM  filter  in  extracting 
the  narrow  band  SF,  spectral  feature  from  the  raw  interferogram 
segment . 


4.2.2  Pattern  Recognition  Considerations. 

The  identification  of  a  specific  IR  vapor  signature  is 
often  made  using  a  pattern  recognition  method  known  as  linear 
discriminant  analysis.  The  linear  discriminant  classifies  a  FT-IR 
interferogram  on  the  basis  of  the  presence  or  absence  of  a 
specific  spectral  signature.  The  linear  discriminant  is  described 
by  equation  2.” 


N 

g(x)  =  w„  +  £  Wi*Xj  (2) 

i=l 

The  Wj  term  denotes  a  weight  vector.  The  response  vector  is 
represented  by  Xj  term.  The  w^  term  is  a  threshold  vector.  The 
g(x)  term  is  the  spectral  response  function.  The  linear 
discriminant,  equation  2,  divides  the  feature  space  with  a 
hyperplanar  surface.  Thus,  the  weight  vector  has  a  number  of 
dimensions  that  are  defined  by  the  interferogram  points. 
Determination  of  the  weight  vector  requires  using  a 
representative  training  data  set.  The  training  set  is  separated 
with  an  optimum  weight  vector  by  using  a  convergence  technique 
based  on  the  perceptron  convergence  theorem.  This  vector  gives  a 
minimum  response  to  spectral  interferants  and  a  maximum  response 
to  the  target  vapor.  The  perceptron  convergence  theorem  uses  an 
iterative  approach  to  determine  the  weight  vector.  It  begins  by 
locating  all  misclassif ications.  These  misclassif ications  are 
used  to  estimate  a  new  weight  vector.  The  process  continues  until 
all  the  data  are  correctly  classified  or  it  is  assumed  that  the 
data  set  is  linearly  inseparable. 

The  constraint  of  linear  separability  is  often  too 
limiting.  The  present  study  uses  a  technique  called  piecewise 
linear  discriminant  analysis  (PLDA).^  This  technique  uses  pieces 
or  portions  of  several  linear  discriminants  (i.e.,  multiple 
linear  discriminant)  to  approximate  a  nonlinear  surface  that 
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better  classifies  the  interferogram  in  pattern  space.  The  PLDA 
approach  also  incorporates  a  simplex  optimization  for  a  more 
effective  placement  of  the  individual  linear  discriminants. 

4.2.3  Digital  Filter  Design. 

Digital  filter  and  pattern  recognition  coefficients 
have  been  generated  for  acetone,  MEK,  and  SF«.  The  most  intense 
spectral  band  for  each  vapor  in  the  atmospheric  window  has 

been  used  in  the  filter  design  (see  Table  2,  column  l) .  A 
description  is  given  below  for  the  three  steps  necessary  to 
generate  robust  sets  of  coefficients.  The  coefficients  are  vapor 
specific.  Description  of  the  procedure  uses  only  MEK,  since  the 
same  technique  applies  to  each  vapor.  Results  of  the  procedures 
that  are  used  for  acetone  and  SF#  are  contained  in  Appendix  C. 

The  first  step  was  to  acquire  a  representative  subset 
from  the  22960  inter ferograms  in  the  full  training  set.  This 
training  subset  contained  10000  inter ferograms  which  included 
1298  vapor  containing  interferograms  and  8702  background 
interferograms.  The  training  subset  selection  was  made  with  the 
global  sampling  algorithm  developed  by  Carpenter  and  Small. A 
tally  of  the  types  of  interferograms  that  were  placed  in  the 
training  subset  is  provided  in  Table  1. 

The  second  step  was  generation  of  the  FIRM  filter 
coefficients.  Smaller  subsets  of  about  400  interferograms  were 
selected  from  the  1298  vapor  containing  interferograms  in  the 
training  subset.  The  FIRM  filter  coefficients  were  generated 
according  to  the  techniques  described  by  Small  and  Harms."  Input 
parameters  important  to  the  performance  of  the  filter  included 
interferogram  segment  length,  interferogram  segment  location,  the 
Gaussian  bandpass  function  peak  position,  and  the  Gaussian 
bandpass  function  FWHH.  The  Gaussian  function  as  noted  previously 
was  used  to  produce  the  ideal  filtered  spectrum.  In  this 
investigation  only  the  filtered  interferogram  segment  length  was 
held  constant  at  75  interferogram  points.  The  75  point  segment 
was  shown  optimal  for  SFj."  The  VOCs  with  wider  spectral  band 
widths  (see  Table  2,  column  2)  might  be  optimal  with  shorter 
interferogram  segment  lengths  (i.e.,  lower  resolution).  Further 
investigation  is  required  to  test  this  hypothesis. 

The  third  step  was  to  determine  the  optimum  for  the 
three  remaining  filter  parameters.  Sets  of  759  different  FIRM 
filters  for  different  Gaussian  peak  positions  were  generated  over 
a  range  of  interferogram  segment  locations  and  band  widths.  To 
assess  the  digital  filter  performance  before  attempting  to 
perform  a  PLDA,  a  one-dimensional  discriminant  was  computed  for 
identifying  potentially  effective  filters.  The  different 
intensity  responses  of  the  filtered  interferogram  segments  in 
Figures  5c  and  5d  suggested  that  the  sum-of -squares  (SOS)  of  the 
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filtered  interferogram  points  would  provide  such  a  one- 
dimensional  discriminant.  Separate  average  magnitude  responses 
were  found  with  the  SOS  calculation  for  subsets  of  the  vapor  and 
background  inter ferograms  in  the  training  set  (see  Table  1) .  The 
ratio  of  vapor  to  background  averaged  responses  is  denoted  as  the 
signal-to-noise  ratio  (S/N) .  A  three  dimensional  plot  of  the  S/N 
for  varying  interferogram  segment  locations  and  Gaussian  function 
band  widths  is  shown  in  Figure  6.  The  axis  for  interferogram 
segment  locations  is  labeled  with  the  starting  point  in  the 
segment  closest  to  the  interferogram  ZPD.  The  maximum  S/N  in 
Figure  6  was  an  interferogram  segment  location  displaced  105 
points  from  the  ZPD  and  a  Gaussian  function  band  width  greater 
than  24  cm'.  The  fourth  dimension,  the  Gaussian  peak  location, 
was  not  shown  in  Figure  6.  The  Gaussian  peak  location  of  1150  cm' 
was  used  in  this  plot.  This  particular  Gaussian  peak  was 
subsequently  found  to  produce  the  lowest  number  of  misclassif ied 
inter ferograms  when  a  PLDA  was  performed.  Thus,  the  filter 
coefficients  were  de  facto  selected  by  the  same  criterion  of 
lowest  misclassif ications.  Multiple  linear  discriminants  were 
generated  with  the  PLDA  approach.'"  The  number  of  inter ferograms 
classified  incorrectly  out  the  training  set  of  10000 
interferograms  is  shown  in  Figure  7  as  a  function  of  Gaussian 
peak  position  and  band  width.  The  lowest  number  of  incorrect 
classifications  was  for  a  peak  position  of  1150  cm'  and  band 
width  of  27  cm'.  These  results  are  for  a  75  point  segment 
displaced  105  points  from  the  interferogram  ZPD.  These  findings 
along  with  those  for  acetone  and  f  are  summarized  in  Table  3. 


Table  3.  Optimal  Filter  Response  Parameters 


INTERFEROGRAM 

POINT 

SEGMENT 

LOCATION 

GAUSSIAN 

PEAK 

POSITION 

(cm') 

GAUSSIAN 

BANDWIDTH 

(cm') 

FILTER  1 

PEAK  1 

RESPONSE  1 
( cm' )  1 

ACETONE 

80  -  155 

1195 

31 

1190  1 

MEK 

105  -  180 

1169 

27 

1150  1 

SF, 

170  -  245 

941 

14 

945  1 

5.  RESULTS  AND  DISCUSSION 

The  performance  of  the  digital  filter  and  multiple 
linear  discriminant  is  best  evaluated  with  the  results  obtained 
from  a  separate  prediction  set  of  interferograms.  Table  1  lists 
the  number  of  interferograms  that  have  been  obtained  with  the 
spectrometer  configurations  of  active  bistatic,  passive 
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Discriminant  Score 


Interferogram  number 


Figure  10.  Discriminant  Scores  for  the  Passive  Laboratory 
Spectrometer  Configuration. 


terrestrial,  and  passive  laboratory  for  the  prediction  test  set. 
The  classification  results  for  the  digital  filter  and  multiple 
linear  discriminant  are  summarized  in  Table  4.  Two  types  of 
misclassif ications  or  incorrect  alarms  occur.  A  false  alarm 
occurs  when  no  vapor  is  present,  but  the  digital  filter  and 
multiple  linear  discriminant  indicates  a  positive  detection.  A 
missed  alarm  results  for  no  detection  when  the  vapor  is  actually 
present.  The  percentages  for  false  and  missed  alarms  are 
determined  from  the  ratios  of  misclassif ications  to  backgrounds 
and  vapor  containing  inter ferograms  respectively  (reference  Table 
1  for  denominator  values  in  ratio) .  The  false  alarms  on  the 
training  set  are  set  to  zero  as  a  computational  constraint.  The 
acetone  and  MEK  results  are  comparable  for  missed  alarms.  The 
slightly  higher  percentage  for  missed  alarms  in  the  prediction 
set  is  to  be  expected.  However,  the  false  alarm  percentage, 
though  reasonable  for  acetone  is  about  an  order  of  magnitude 
higher  for  MEK  and  comparable  in  magnitude  to  the  missed  alarms. 
If  vapor  containing  inter ferograms  are  incorrectly  classified  in 
the  training  set  then  a  higher  false  alarm  count  results.  A 
closer  examination  of  the  training  set  is  required  to  test  this 
possibility. 

The  classification  results  for  SF^  are  also  listed  in 
Table  4.  The  number  of  interferograms  are  listed  in  Table  1  that 
are  selected  for  training  and  prediction  subsets.  These  subsets 
differ  in  that  the  prediction  subset  contains  only  passive 
laboratory  data  for  the  vapor-present  case.  On  the  other  hand, 
the  SFj  training  subset  includes  interferograms  from  all  three 
spectrometer  configurations  for  the  vapor-present  case.  Using 
only  the  passive  laboratory  data  for  the  prediction  subset 
provides  a  more  accurate  determination  of  the  causes  for  false  or 
missed  alarms.  The  number  of  false  alarms  for  SF«  is  comparable 
to  NEK.  However,  closer  examination  of  the  data  reveals  about  20 
of  the  25  false  alarms  are  made  for  conditions  near  or  less  than 
the  detection  limit  of  the  spectrometer.  Using  a  reference 
blackbody  background  subtraction  indicates  that  only  about  5  of 
the  20  false  alarms  are  valid.  Therefore,  a  more  realistic  count 
of  false  alarms  is  10  (i.e.,  0.5%). 

The  FIRM  filter  parameters  are  given  in  Table  3.  These 
parameters  have  been  used  in  the  generation  of  the  Table  4 
classification  results.  A  plot  of  the  FIRM  filter  frequency 
responses  is  shown  in  Figure  8.  The  frequency  response  plots  are 
generated  with  the  approach  that  is  given  in  Appendix  D.  The 
attenuation  provided  by  the  FIRM  filters  is  consistent  with  the 
signal  dynamic  range.  For  interferogram  segments  displaced  off 
the  interferogram  center  burst  only  about  nine  bits  of  the 
analog-to-digital  converter  range  is  used.  Therefore,  a  filter 
attenuation  in  excess  of  about  27  decibels  (dB)  is  unnecessary. 
The  attenuation  in  the  filter  bandpass  removes  the  residual 
broad  band  detector  envelope.  Hence,  the  bandpass  attenuation 
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becomes  greater  as  the  Interferogram  segment  location  approaches 
the  interferogram  center  burst  region.  The  gradual  decrease  of 
attenuation  away  from  the  filter  bandpass  reflects  the  inherent 
optical  filtering  of  the  detector  response  curve  (i.e.,  detector 
cutoffs  near  800  cm'  and  1600  cm') . 


Table  4.  Classification  Results 


FALSE 

MISSED 

TOTAL  1 

ALARMS 

ALARMS 

CORRECT 

Acetone 

0 

27 

9973 

(0.0%) 

(2.6%) 

(99.7%) 

TRAINING 

RESULTS 

MEK 

0 

(0.0%) 

63 

(4.8%) 

9937 

(99.4%) 

SFs 

0 

11 

9989 

(0.0%) 

(1  2%) 

(99.9%) 

Acetone 

16 

268 

9716 

(0.3%) 

(5.5%) 

(97.2%) 

PREDICTION 

RESULTS 

MEK 

259 

(3.9%) 

201 

(6.0%) 

9540 

(95.4%) 

SF* 

25 

89 

9886 

(1.3%) 

(1.1%) 

(98.9%) 

The  optimal  FIRM  filter  and  multiple  linear 
discriminant  performance  occurs  for  the  minimum  in  classification 
errors  for  the  training  subset.  The  optimum  for  the  VOCs  gives  a 
minimum  error  classification  count  for  a  filter  response  curve 
centered  off  the  spectral  band.  The  peak  filter  bandpass 
responses  differ  from  the  acetone  and  MEK  spectral  bands  by 
26  cm'  and  22  cm'  respectively.  No  such  difference  is  discernible 
for  the  SF«  filter  response.  The  SF^  case,  unlike  the  VOCs  does 
not  possess  a  lower  intensity  band  adjacent  to  the  major  spectral 
feature.  This  lower  intensity  band  becomes  important  for  the 
bistatic  active  spectrometer  configuration.  Figure  9  shows  the 
absorbance  spectra  for  acetone  and  MEK  that  is  taken  in  the 
active  bistatic  spectrometer  configuration.  The  peak  spectral 
positions  are  delineated  with  solid  lines,  while  the  broken  lines 
indicate  the  peak  filter  pass  band  response.  It  is  important  to 
note  that  the  digital  filtering  of  the  short  interferogram 
segment  is  only  one  step  in  the  background  suppression  method. 

The  classification  is  made  with  the  multiple  linear  discriminant 
that  is  generated  by  PLDA^  on  the  filtered  interferogram  segment. 
The  multiple  linear  discriminant  uses  pertinent  band  contour 
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Information  that  facilitates  a  separation  in  the  pattern  space. 

In  the  VOC  case,  frequency  components  from  the  lower  intensity 
band  that  is  adjacent  to  the  major  spectral  band  become 
important.  These  results  indicate  that  more  work  is  needed  to 
assess  the  importance  of  functional  filter  forms  other  than  the 
Gaussian  for  direct  interferogram  processing. 

The  VOC  FIRM  filter  and  multiple  linear  discriminants 
have  been  implemented  on  an  IBM/PC  compatible  computer.*^  The 
results  of  the  real-time  operation  with  the  MIDAC  spectrometer  is 
shown  in  Figure  10.  Figure  10  is  a  plot  of  the  multiple  linear 
discriminant  score  as  a  function  of  sequentially  collected 
inter ferograms.  These  interferograms  are  collected  in  the  passive 
laboratory  spectrometer  configuration.  A  fixed  concentration  of 
MEK  is  contained  in  the  vapor  cell.  The  blackbody  background 
temperature  is  varied  around  the  cell  ambient  temperature  by 
approximately  5*’C  increments  for  a  blank  cell,  vapor  cell,  and  no 
cell.  The  distribution  of  interferograms  are  as  follows:  0-97, 
blank  cell  background;  98-187,  vapor  cell;  and  188-278,  no  cell. 
The  background  and  vapor  cell  temperatures  are  equal  for 
interferograms  141  to  155.  The  discriminant  performs  as  expected 
with  no  response  for  interferograms  141  to  155.  Clearly,  a 
nonzero  temperature  differential  between  the  vapor  and  the 
background  is  required  to  produce  a  spectral  emission  or 
absorption  vapor  spectrum. 

6 .  CONCLUSIONS 

The  application  of  FIRM  filters  on  selected 
interferogram  segments  followed  by  a  multiple  linear  discriminant 
analysis  provides  a  means  of  suppressing  the  broad  band  detector 
envelope.  This  removes  the  need  for  a  stable  reference  background 
spectrum.  Often  a  background  reference  spectrum  does  not  exist 
for  many  open-air  monitoring  conditions,  especially  with  passive 
terrestrial  spectrometer  configurations.  The  background 
suppression  approach  is  found  to  work  well  in  the  discrimination 
for  the  vapors  studied,  while  rejecting  spectral  information  from 
other  atmospheric  constituents.  This  discrimination  capability  is 
essential  for  many  atmospheric  monitoring  applications. 

Acetone,  MEK,  and  SF^  provide  spectral  bands  of 
decreasing  bandwidth.  These  bandwidths  are  shown  to  relate 
directly  to  the  optimxim  interferogram  segment  selected  for 
digital  filtering  with  the  broader  bands  residing  nearer  the 
center  burst.  The  multiple  linear  discriminant  analysis  has 
clearly  demonstrated  the  importance  of  band  contours  when 
contrasting  the  results  for  the  VOCs  versus  SF^.  Direct 
processing  of  interferograms  with  digital  filtering  combined  with 
multiple  linear  discrimination  facilitates  the  use  of  FT-IR 
spectrometry  in  open-air  monitoring  applications. 
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APPENDIX  A 


DATA  COLLECTION  DISK  FORMAT 


GLOBAL  HEADER  >  512  BYTES 


Byte 

Bytes 

Field 

Data 

Position 

Used 

Description 

Type 

«General  Information» 

0 

10 

Filename  (CCC/###) 

String 

10 

10 

Date  (MM/DD/YY) 

String 

20 

10 

Start  Time  (HH:MM:SS) 

String 

30 

10 

Stop  Time  (HH:MM:SS) 

String 

40 

2 

Stop  Scan  Number 

Integer 

42 

10 

Operator's  Name 

String 

52 

44 

Unused 

«Sensor  Information» 

String 

96 

20 

Sensor  Identification 

String 

116 

2 

Collection  Mode 

Integer 

118 

2 

Integer  Type 

Integer 

120 

2 

Interferogram  Pts/Scan 

Integer 

122 

8 

Resolution 

Double 

130 

8 

Scan  Speed 

Double 

138 

8 

Mirror  Velocity 

Double 

146 

8 

Sampling  Frequency 

Double 

154 

8 

Starting  Frequency 

Double 

162 

8 

Ending  Frequency 

Double 

170 

8 

Maximum  Wavenumber 

Double 

178 

2 

HeNe  Zero  Crossings/Pt. 

Integer 

180 

16 

Unused 

«Weather  Information» 

String 

196 

2 

Ambient  Temperature 

Integer 

198 

8 

Barometric  Pressure 

Double 

206 

2 

Humidity 

Integer 

208 

2 

Wind  Speed 

Integer 

210 

2 

Wind  Direction 

Integer 

212 

2 

Sensor  FOV  Direction 

Integer 

214 

2 

Precipitation  Code 

Integer 

216 

40 

Unused 

«Comments» 

String 

256 

64 

comment  -  line  #1 

String 

320 

64 

comment  -  line  #2 

String 

384 

64 

comment  -  line  #3 

String 

448 

64 

comment  -  line  #4 

String 
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SUBFILE  INTERFEROGRAM  HEADER  -  64  BYTES 


Byte 

Bytes 

Field 

Data 

Position 

Used 

Description 

Type 

0 

2 

Interferogram  number 

Integer 

2 

10 

Filename 

String 

12 

10 

Time  (HH:MM:SS) 

String 

22 

2 

Center  burst  location 

Integer 

24 

2 

ADC  Gain 

Integer 

26 

2 

Number  of  Averaged  Intf. 

Integer 

28 

34 

Unused 

String 

62 

2 

Error  Code 

Integer 

INTERFEROGRAM  POINTS  -  VARIABLE  BYTES 


Integer  data  block  consists  of  consecutive  data  bytes  with  data 
type  designated  in  the  global  header.  The  integer  data  block 
length  is  specified  by  the  number  of  interferogram  points  per 
scan  (see  global  header) . 
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APPENDIX  B 


VAPOR  TRAINING  SUBSET  DESCRIPTION 
TRAINING  SET  DATA  -  Acetone  Vapor  Present  - 


ACTIVE  BISTATIC  (850) 


IR  Source 

* 

Location 

Weather 

FOV 

1.0  V 

112 

Site  3 

Rain 

NT 

1.5  V 

123 

Site  2 

Fair 

T 

2.0  V 

147 

Site  3 

Rain 

NT 

3.0  V 

121 

Site  3 

Rain 

NT 

4.0  V 

110 

Site  3 

Rain 

NT 

5.0  V 

63 

Site  3 

Rain 

NT 

6.0  V 

64 

Site  3 

Rain 

NT 

8.0  V 

48 

Site  2 

Fair 

T 

12.0  V 

62 

Site  2 

Fair 

T 

PASSIVE  TERRESTRIAL  (89) 


IR  Source 

Location 

Weather 

FOV 

Shed 

Site  3 

Fair 

NT 

Shed 

8 

Site  3 

Rain 

NT 

Shed 

24 

Site  3 

Rain 

NT 

Shed 

48 

Site  3 

Rain 

NT 

PASSIVE  LABORATORY  (84) 

Concentration 

aT 

# 

Vapor 

(»C) 

14 

XI 

29 

16 

XI 

19 

15 

XI 

9 

22 

XI 

-1 

17 

XI 

-11 

Note:  Superscript  on  the  inter ferogram  count  number  (#) 

indicates  the  number  of  inter ferograms  averaged 
per  interferogram  subfile. 
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TRAINING  SET  DATA  -  MEK  Vapor  Present 


ACTIVE  BISTATIC  (1121) 


IR  Source 

II 

11 

II 

•tfe  It 
II 
II 
II 
II 
II 
II 

Location 

Weather 

FOV 

1.0  V 

202 

site  3 

Cloudy 

NT 

1.0  V 

75 

Site  3 

Rain 

NT 

1.5  V 

56 

Site  4 

Fair 

T 

2.0  V 

94 

Site  3 

Cloudy 

NT 

2.0  V 

199 

Site  3 

Rain 

NT 

3.0  V 

142 

Site  3 

Rain 

NT 

4.0  V 

136 

Site  3 

Rain 

NT 

5.0  V 

71 

Site  3 

Cloudy 

NT 

5.0  V 

51 

Site  3 

Rain 

NT 

6.0  V 

71 

Site  3 

Cloudy 

NT 

6.0  V 

24 

Site  3 

Rain 

NT 

PASSIVE  TERRESTRIAL 

(103) 

IR  Source 

Location 

Weather 

FOV 

Shed 

20 

Site  3 

Fair 

NT 

Shed 

26 

Site  3 

Cloudy 

NT 

Shed 

4« 

Site  3 

Cloudy 

NT 

Shed 

53 

Site  3 

Rain 

NT 

PASSIVE  LABORATORY 

(74) 

i 

Concentration 

Vapor 

aT 

(»C) 

15 

XI 

29 

15 

XI 

19 

12 

XI 

9 

16 

XI 

-1 

16 

XI 

-11 
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TRAINING  SET  DATA 


-  SFj  Gas  Present 


ACTIVE  BISTATIC  (211) 


IR  Source 

1 

1 

1 

1 

Location 

Weather 

FOV 

8.0  V 

120 

Site  1 

Rain 

T 

12.0  V 

91 

Site  2 

Fair 

T 

PASSIVE  TERRESTRIAL  (57) 

IR  Source 

* 

Location 

Weather 

FOV 

Shed 

33 

Site  3 

Fair 

NT 

Shed 

24'* 

Site  3 

Fair 

NT 

PASSIVE  LABORATORY  (660) 


Concentration 

aT 

1 

1 

1 

1 

(cc) 

(“C) 

15 

0,2 

29 

16 

0.2 

19 

13 

0.2 

9 

16 

0.2 

4 

23 

0.2 

1 

16 

0.2 

-11 

100 

0.04 

29 

100 

0.04 

19 

61 

0.04 

9 

26 

0.04 

-1 

100 

0.04 

-11 

100 

0.08 

7 

60 

0.08 

5 

14 

0.08 
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APPENDIX  C 


ACETONE  AND  SFj  PLDA  SUMMARY 


ACETONE  FILTER  S/N  (1187.7  WAVENUMBERS) 


ACETONE  PATTERN  RESULTS 
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APPENDIX  D 


FIRM  FILTER  FREQUENCY  RESPONSE 


1.  Select  a  spectral  frequency,  between  0  and 

fj  «  2»v*Pj  -  interferometer  modulation  frequency 
V  »  interferometer  mirror  velocity 

2.  Input  the  synthetic  cosine  waveform  for  the  selected  spectral 

frequency. 

Xjj  -  cos(27ri(fj/f^) ) 

Xjj  »  simulated  interferogram  point  at  frequency  fj 
i  *  interferogram  point  index 
fj  *  selected  input  frequency 

f^  s  highest  frequency  response  to  be  calculated 

3.  Calculate  filter  response  on  the  logarithmic  decibel  scale  . 

(dB)  using  SOS  of  filtered  and  unfilter  synthetic 
interferogram  segments 

dB  *  -20»log(yj/Xj) 

N 

Xj  =  E  input  unfiltered  segment  SOS 

i=l 

N 

Yj  =  E  (Yy)^  output  filtered  segment  SOS 
i=l 

4.  Oscillatory  behavior  of  frequency  response  plots  due  to  the 

discrete  number  of  frequency  components  used  to  generate  the 
frequency  response  and  discontinuous  interferogram 
segments  used  by  the  FIRM  filter. 
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